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the conflict matrix for this investigation. 

Specific investigation of the model's behavior is observed 
varying the Command, Control, Communications and Intelligence 
(C3i) enhancements of one combat force against that of an 
opposing force while one force utilizes misinformation and 
deception as a counter C3 tactic and the other uses physical 
destruction as a counter C3 tactic. 

Potential military use of this model as an analytical tool 
for playing "what if?" type wargames is envisioned following 
further research and study of the effects on model behavior by 
varying other parameters not specifically addressed in this 
initial investigation. 
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ABSTRACT 

The Dynamic Model of Modern Military Conflict developed 
by Dr. Paul Moose, Naval Postgraduate School, is described 
by its system of differential equations followed by an inves- 
tigation of its behavior. This investigation is predicated 
by an analysis of the model's stability about equilibrium 
using a method attributed to the study of ecosystems. The 
basis for this analysis is the formulation and subsequent 
evaluation of a community matrix termed the conflict matrix 
for this investigation. 

Specific investigation of the model's behavior is observed 
varying the Command, Control, Communications and Intelligence 
(C^I) enhancements of one combat force against that of an 
opposing force while one force utilizes misinformation and 
deception as a counter C^ tactic and the other uses physical 
destruction as a counter C^ tactic. 

Potential military use of this model as an analytical 
tool for playing "what if?" type wargames is envisioned 
following further research and study of the effects on model 
behavior by varying other parameters not specifically addressed 
in this initial investigation. 
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INTRODUCTION 



I . 

A military conflict has a variety of elements that 
potentially influence its eventual outcome. Commanders of 
the various echelons within the combat environment determined 
by a particular military conflict are faced with decisions 
that require the fusion of information from countless sources. 
These sources may be highly sophisticated sensors, relaying 
timely digital information concerning enemy force concentra- 
tions or movements, aircraft or missile tracks both friendly 
and nonfriendly, or perhaps intercepted communications signals 
These sources could also be the professional schools, exercise 
participation or intelligence documents that have molded the 
commander's decision process over his military career whether 
it be extensive or minimal. Regardless of the case, decisions 
are made with some reliance on or rejection of prior knowledge 
The assimulation of this prior knowledge, whether relative to 
the situation or not can be a monumental task in itself. 

The sophistication of computer driven information systems 
and a presumed need for rapid decision making either by the 
man as a commander or the machine as his extension has caused 
a proliferation of model development in recent years [Ref. 1] . 
These models attempt to reflect the variety of elements that 
potentially influence the outcomes of military engagements, 
manipulate them according to some predetermined scheme or 
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algorithm and finally output some prescription or calculated 
result. The output of such, models, together with the parti- 
cular input elements and scheme of manipulation offer 
limitless applications either for direct command decision 
making, training or running countless "what if" exercises 
relative to known or suspected enemy and friendly combat 
information or characteristics. 

A. PURPOSE 

The primary purpose of this thesis is to initiate inves- 
tigation of one such model that exhibits the potential for 
use in a variety of areas, but in particular in looking at 

the "what if" questions associated with command, control and 

-5 3 

communications (C-^) and counter C . The specific model under 
investigation is the Dynamic Model of Modern Military Conflict 
proposed by Dr. Paul Moose, Academic Associate for the Joint 
Command, Control and Communications curricula at the Naval 
Postgraduate School [Ref. 2]. 

B. SCOPE 

The basic description of the model is given in Chapter 2. 
This description is in the form of a system of differential 
equations similar to that which describes Lanchester-type 
attrition models, however it includes terms relative to C^ 

"3 

and counter C for opposing forces. In this respect, informa- 
tion possessed by opposing forces is elevated to an equivalent 
level with the forces themselves. Furthermore, the description 
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includes a discussion on how the equations of this model are 
also similar to evolution equations which, describe multi- 
species ecosystems. Such, systems include natural death (loss) 
rates and replenishment of species within a population thus 
leading to a mathematical discussion of this model's behavior 
relative to an established equilibrium. 

Following this description. Chapter 3 looks at the 
stability of the model in terms of a "community matrix", 
another reference to ecosystems which incorporates those 
parameters which tend to increase or reduce species within 
a population or community. This "community matrix" is 
renamed the "conflict matrix" for this model and explored 
by assigning a certain set of asymmetrical values to the 
elements of the matrix relative to hypothetical opposing 
forces in a military conflict. Specifically the investiga- 
tion undertaken in this chapter, demonstrates conditions 
where opposing forces stabilize at a set equilibrium point 
or diverge from this point in unstable fashion. 

Chapter 4 looks at more specific behaviors relative to 
the stability and instability demonstrated in Chapter 3. 

The identical initial conditions used in Chapter 3 are 
again used in this chapter and represent arbitrary units 
rather than specific calculated values. These values were 
chosen purely to demonstrate the general behavior of the 
model rather than replicate any specific military conflict. 
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Finally effects on equilibrium relative to the initial 
conditions used during stability and behavior investigations 
are explored in Chapter 5. Specifically, the initial condi- 
tions are aimed at isolating the effects of varying the 
command, control, communications and intelligence (C^I) 

3 

enhancements of one force against an opposing force's C I 
enhancements. Utilizing these results, further demonstration 
of the model behavior is then achieved by varying resource 
and information variables associated with opposing military 
forces . 

In essence, this thesis is concerned with an initial 

demonstration that the Dynamic Model of Modern Military 

Conflict can perform as might be expected rather than prove 

that it works in all cases. This approach opens several 

avenues to continued research leading to a more complete 

demonstration of its usefulness in evaluating the influence 
3 3 

of C I and counter C in modern military conflict. 
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II . A DYNAMIC MODEL OF MODERN MILITARY CONFLICT 
A. DEFINITION 

The Dynamic Model of Modern Military Conflict incorporates 



those found in Lanchester-type models. Specifically the 
effects of the above are described in terms of information 



opposing X and Y combat forces [Ref. 3]. The system of 
differential equations which characterize this model are 
given below followed by sections in this chapter which 
briefly touch on the background of the model and provide 
specific definitions of the terms involved. 



3 

the effects of Cl with basic attrition coefficients such as 



(I and I ) and general purpose resources (M and M ) for 
X y ^ Y 



(Eqn 2.1) 







(Eqn 2.2) 




(Eqn 2.3) 



d(I )/dt = I [-A - a I - Y M ] + C I M + 

y y y yx X 'yx X-* yy ye y 



yy ye y 




(Eqn 2.4) 



d(M,)/dt = M [-3 - 6 I - B M ] - d 

y y y yx x yx x-* 
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B . BACKGROUND 



This system of four first order non-linear differential 

d ( S ) 

equations is of the form - = _p(S) + ^ set of 

equations in this form are known as evolution equations. 

S and Q are 4x1 column vectors and there are four func- 
tional relationships which are considered to be at most 
quadratic. Expanding the loss function term gives the 
following ; 





(Eqn 2.5) 
4 

^ a..S./ i — 1,2, 3/4 

j=i : 



i .e . 



- (S)= - quadrative term - linear term (Eqn 2.6) 

This system as it turns out is a very generalized form of the 
Lotka-Vol terra equations for multi-species ecosystems [Ref. 

4: 37-38]. More specifically it takes on the form: 

(Eqn 2.7) 



d(S^) 

dt 



-s. [ I 
" jil 



a. .S.j- i 
13 3 > 



^ a. .S . + Q. , 
j=l 3 



i = 1,2, 3, 4 



C. EXPLANATION OF TERMS 

In general the greek letter coefficients are quadratic, 
and the others are linear. Following is the definition of 
the terms that characterize the model: 

I - Information of force X 
I„_ - (Same as above at equilibrium) 

M - General purpose resources for force X 
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M 



xe 



X 



B 



yx 



yx 



X 



X 



XX 



XX 



a 



yx 



Y 



yx 



yx 



yx 



- (Same as above at equilibrium) 

3 

- Natural death (loss) rate for C I resources 
of force X 

- Natural death (.loss) rate for general purpose 
resources of force X 

- Lanchester "area fire" coefficient, X on Y 

- Lanchester "aimed fire" coefficient, X on Y 

3 

- Fixed replenishment rate for C I resources of 
force X 

- Fixed replenishment rate for general purpose 
resources of force X 

- Replenishment from general purpose resources 
for C3i production, force X 

3 

- Diminishment of X force for C I production 

3 

- Counter C by force X on force Y by misinformation 
and deception from C3i resources 

- Counter C^ by force X on force Y by physical 
destruction from general purpose resources 

3 

-Cl (quadratic) force effectiveness enhancements, 

X on Y 

- C^I (linear) force effectiveness enhancements, 

X on Y 



The remaining coefficients with y subscripts or xy subscripts 
refer to the Y force or Y on X activities. 



D. EXPECTED MODEL BEHAVIOR 

The actual behavior of the model is dependent initially 

on the equilibri\im for the system. One equilibrium point is 

determined by solving the following linear set of equations 

for I ^ M , I and M : 

-^xe' xe ye ye 
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(Eqn 2.8) 
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Once equilibrium is established the effects of varying 

certain coefficients on dependent variables I , M , I , M 

x' x' Y Y 

can be illustrated by solving the system of differential 
equations over the independent variable time. 
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III. STABILITY ANALYSIS 



A. DEFINITION 

Stability, as it relates to a deterministic model or 
deterministic system of equations, is more specifically 
defined as "neighborhood stability" or "stability in the 
vicinity of equilibirum" . For population models with 
environmental parameters all well defined constants, 
interest centers on the community equilibrium where all 
species' populations have time-independent values, i.e. 
all net growth (death) rates are zero. Such an equilibrium 
may be called stable if, when the populations are perturbed, 
they in time return to their equilibrium values; the return 
may be achieved either as damped oscillations or monotonically . 
Conversely, if such a disturbance tends to amplify itself, the 
system is called unstable; such instability may appear as an 
oscillatory or as a monotonic growth in the disturbance. 

This definition is sufficient for linearized systems, however, 
it may be misleading for non-linear systems when full global 
stability of a system is of concern [Ref. 4:15]. 

Full global stability is said to be characterized by 
neighborhood analysis if the existence of a Lyapunov function 
can be determined [Ref. 4:15]. The existence of such a 
function is not explored in this thesis, however, regardless 
of whether such a function exists or not, neighborhood analysis 
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is within the scope of this thesis and founded in the 
generation of a community matrix. 



B. COMMUNITY MATRIX 

The community matrix is said to both summarize the system 
(its elements being determined by the interactions between and 
within species near equilibrium) and sets the neighborhood 
stability by the sign of its eigenvalues (all negative real 
parts for all eigenvalues imply stability) . Specifically it 
is generated by the following mathematical derivation [Ref. 4; 
19-22] . 

Since the model is of the form 



d(S) 

~dF" 



= -F(S) + Q 



its multispecies population dynamics are given by the set 
of four equations; 

(Eqn 3.1) 



d(S^(t) ) 
dt 



= F^(S^(t), S2(t), S3(t), S^(t)), i = 1,2,3, 



■The growth rate of the i^ species at time t is given by 

the nonlinear function of all relevant interacting 

* 

populations. Equilibrium populations, S^ , follow from 
the 4 algebraic equations obtained by setting all growth 
rates at zero: 

(Eqn 3.2) 

★ ★ -k ★ 

0 = F^(S^ , S2 , S3 , S^ ) 
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(Eqn 3.3). 



For each, population, let 



S^(t). = + s^(t)_ 

where measures the initially small perturbations to 

the i^ population. Expanding each of the basic equations 
around this equilibrium in a Taylor series and discarding all 
second order and higher terms in the population perturbations 
s, a linearized approximation is obtained: 



d s^(t) 
~Ht 



= y a . . S . ( t) 



(Eqn 3.4) 



j = l 



Equivalently, in matrix notation 

= A X(t). 

where X is the 4x1 column Matrix of the 



’i' 



(Eqn 3.5) 
and A is 



the four by four "community matrix" whose elements a^^ 
describe the effect of species j upon species i near 
equilibrium. The elements a^^ depend both on the details 
of the original equations and on the values of the equilibrium 
population through 



3F. * 

a . . = - , 1, 
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^3S 

3 



(Eqn 3.6) 



The partial derivatives, denoting the derivatives of F^ 
keeping all populations except constant, are then 

evaluated with all populations at their equilibrium values. 
The corresponding community matrix for the Dynamic Model of 
Modern Military Conflict is called the "conflict matrix" 
and its elements are given in the following matrix: 
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-(A + a + Y ) 
• X xy 



XX 



-a 



yx 

-(6 +d ) 
yx yx 



XX 



-a 



xy 



-Y. 



xy 



-(B +6 +3 i -(6 +d ) -C3 +b ) 

X xy xy xy xy xy xy 



~Y 

'yx 

-(3 +b ) 
yx yx 



-(A +a +Y ) c 
' y yx 'yx' yy 



yy 



-(B +5 +3 ) 

y yx yx 



C. METHOD OF ANALYSIS 



To investigate the stability of the model, eigenvalues 

are computed relative to the conflict matrix given above and 

a set of initial values for the coefficients that comprise 

the elements of the matrix. Two cases are investigated: 

one in which element d is held constant and element d 

yx xy 

is varied across a range of values and another where d^^ 
is held constant and d^^ is varied across a range of values. 
The particular method of computing eigenvalues is the IMSL 
routine EIGRF [Ref. 5j . 



D. INITIAL CONDITIONS 

The initial conditions for investigating the two cases 
are given in Table I. These initial conditions were chosen 
to illustrate some of the potential behavior that this model 

3 

may exhibit. The variable parameters are C I enhancements 
for both X and Y forces. Furthermore the model is 
asymmetrical; X is using physical destruction by forces as 
a counter technique and Y is using misinformation and 

deception as a counter technique. 
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TABLE I 



INITIAL CONDITIONS FOR STABILITY INVESTIGATION 



A=A=B=B= Q.5 

X y X y 

c = c =d =d =0.0 

XX yy XX ^yy 

5 =3 = 1.0 

xy yx 

6 =3 = O.Q 

xy yx 

b = b =0.0 

xy yx 

= 1.00 (misinformation and deception 

V = 

^ xy “ ® ^ 

= 1.0 (physical destruction) 

3 

d = variable C I enhancements (Y on X) 
xy 

3 

d = variable C I enhancements (X on Y) 
yx 

Equilibrium is taken to be the vector 1 , i.e 



I =M =I =M =1 
xe xe ye ye 
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Substitution of the initial parameter conditions into 



the conflict matrix gives the following matrix: 



-1.5 

0 

0 

-d 

yx 



Q 

-1.5 

-l.Q 

-i.a 



- 1.0 

-d 

yx 

-1.5 

0 



0 

- 1.0 

0 

-1.5 



The particular cases investigated then are: 

1 . 

by 0.05. 

2 . 

by 0.05. 



d held constant at 1 and d varied from 0 to 2 
yx xy 



d held constant at 1 and d varied from 0 to 2 
xy yx 



E. RESULTS AND DISCUSSION 

The eigenvalues generated by case 1 are presented in 

Table II. It can be seen that for values of d greater 

xy 

than 0.80 a positive real eigenvalue appears. This is 
indicative of instability. Similarly all real parts of 
eigenvalues for d^^ less than or equal 0.80 are negative 
thus indicating stability. 

The eigenvalues generated by case 2 are presented in 
Table III. For values of d^^ greater than 0.55 instability 
is indicated while stability is indicated for values of d 
less than or equal to 0.55. 
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EIGENVALUES FOR VARIABLE C^l ENHANCEMENTS (Y ON X) 
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EIGENVALUES FOR VARIABLE C l ENHANCEMENTS (X ON Y) 
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The results of the above two cases are summarized 
graphically in Figures 1, 2, and 3. Figure 1 depicts the 
generation of eigenvalues with, positive real parts 
(instability) beyond 0.5 and 0.8 for d^^ and d^^ 
respectively, and negative real parts elsewhere (stability) . 
Figure 2 represents the other real eigenvalue which is 
always negative thus not contributing to instability. Figure 
3 represents the positive imaginary part of the remaining two 
eigenvalues which turn out to be complex conjugate pairs with 
real part equal to -1.5. Existence of these imaginary parts 
indicate that the system does in fact oscillate while dampen- 
ing out. 

The eigenvalues generated by these two cases lead to 
more specific investigation of stability and instability. 

The eigenvalues have located the regions of stability and 
instability so now by using the appropriate values of d 

xy 

and d the behavior of the model can be observed relative 
yx 

to perturbations of the dependent variables over time for 
stable and unstable cases. 
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Real 




1. Eigenvalues with, positive real parts 
(imaginary part =*Q) . 



27 



Real 




Figure 2. Eigenvalues with all negative real 
parts (imaginary part = 0) . 
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Imaginary 

0 12 3 4 




Figure 3. Positive imaginary parts of complex 
conjugate pairs (real part = -1.5). 
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IV. SPECIFIC STABILITY ANALYSIS 



A. METHOD OF ANALYSIS 



More specific illustrations of stability in the Dynamic 
Model of Modern Military Conflict are achieved by perturbing 
the dependent variables and observing the behavior of the 
model relative to time. The initial values used in the 
stability analysis of Chapter 3 are substituted into the 
differential equations which characterize the model. These 
differential equations turn out to be the following for the 
parameter values of Table I. 



d I 



X 



dt 
d M 



X 



dt 
d I 

i 

dt 

d M 

i 

dt 



= 1.5 - I^ (I^ + 0.5) 



= 1.5-d I -M (M + 0.5) + d 

xy y X y xy 



= 1.5 - I^ (M,^ + 0.5) 



= 1.5 - d I - M (M + 0.5) + d 

yx X y X yx 



(Eqn 4.1) 
(Eqn 4.2) 
(Eqn 4.3) 
(Eqn 4.4) 



-The Interactive Ordinary Differential Equations package 

(lODE) was employed using these differential equations and 

representative values for d and d (stable and 
^ xy yx 

unstable values) [Ref. 6]. Specifically, plots were generated 
to graphically represent two behaviors: 

1. Dependent variables against the independent variable 

i.e.,I ,M,I ,M against time. 

X X Y Y 

2. Phase plots of opposing forces i.e., M^ against M^. 
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B. STABILITY WITH AT a. 4 

Table IV shows the initial values for the dependent 

variables, time duration and increment for plotting for 

each of fourteen trial runs of this model. Both d and 

xy 

d^^ were held constant at Q.4 and 1.0 respectively. In 

addition, the final column in Table IV represents the time 

at which all dependent variables damped out to or very nearly 

to the established equilibrium of 1. 

Graphically, Figure 4 depicts the general plot where 

only one dependent variable was perturbed (trials 1 through 

4) . Figures 5 through 8 further illustrate phase plots for 

these same trials of M vs M illustrating return to 

X y ^ 

equilibrium over time. 

Two other trials are illustrated in Figures 9 through 12. 

Specifically Figures 9 and 10 represent perturbations of both 

I and I (trial 6) . The interesting phase plot in Figure 10 
X y 

is indicative of the oscillatory nature of both and 

relative to these perturbations. Figure 11 depicts a triple 
perturbation involving I^, I^ and M^, with the resulting 
oscillatory nature of and shown in Figure 12 (Trial 13) . 

Although the above illustrations are not all inclusive, it 
should be noted that all dependent variables did in fact damp 
out over the set time period and phase plots of and 

revealed return to equilibrium in all trials. 
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TABLE IV 



INITIAL CONDITIONS WITH 







c 


1 

xy 


AT 0.4 AND 


^yx 


1.0 


(STABLE) 






Trial # 


Time 


Increment 


5c 


M 

X 


I 


M 

_y. 


Time for 
9 0 % Damp 
out (sec) 


1 


0 


to 


50 


.05 


1.25 


1.00 


1.00 


1.00 


7 


2 


0 


to 


50 


.05 


1.00 


1.25 


1.00 


1.00 


12 


3 


Q 


to 


50 


.05 


1.00 


1.00 


1.25 


1.00 


9 


4 


Q 


to 


50 


.05 


l.QQ 


1.00 


1.00 


1.25 


10 


5 


0 


to 


50 


.05 


1.25 


1.25 


1.00 


1.00 


17 


6 


0 


to 


50 


.05 


1.25 


1.00 


1.25 


1.00 


1 


7 


0 


to 


50 


.05 


1.25 


1.00 


1.00 


1.25 


1.5 


8 


0 


to 


50 


.05 


1.00 


1.25 


1.00 


1.25 


1.5 


9 


0 


to 


50 


.05 


1.00 


1.00 


1.25 


1.25 


15.5 


10 


0 


to 


50 


.05 


1.00 


1.25 


1.25 


1.00 


4.0 


11 


0 


to 


50 


.05 


1.25 


1.25 


1.25 


1.00 


12.0 


12 


0 


to 


50 


.05 


1.25 


1.25 


1.00 


1.25 


11.0 


13 


0 


to 


50 


.05 


1.25 


1.00 


1.25 


1.25 


12. Q 


14 


0 


to 


50 


.05 


1.00 


1.25 


1.25 


1.25 


4,0 
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Figure 4. Single Perturbation, stable (general plot) . 



t ^66 




fO 

r* 

o 



wn 

rr* 



34 



Figure 5. Single Perturbation, stable phase plot (I^ at 1.25) 
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Figure 6. Single Perturbation, stable phase plot (M^ at 
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Figure 7. Single Perturbation, stable phase plot (I^ at 1.25). 
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Figure 8. Single Perturbation, stable phase plot (M at 
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Figure 10. Double Perturbation, stable phase plot (I^ and 
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Figure 12. Triple Perturbation, stable phase plot (I^, I and 



C. STABILITY WITH AT Q.ia 

Similar to Table IV, Table V lists the initial parameters 

for investigation of the^ model with d at 0.10 and d 

yx xy 

at 1.00. Graphically, for a single perturbation, this form 
of the model produces similar results as before. The excep- 
tion is given in Figure 13 where, only I^ is perturbed 
(trial 3) . and can be seen to return to equilibrium 

without oscillation as had been indicated in the previous 
investigation. Again all trials exhibit the dependent 
variables damping out in the time allowed and since no 
significantly differently graphical forms were produced they 
are not included in this section. 

D. INSTABILITY WITH d^^ AT 1.4Q 

Table VI lists the initial parameters for the unstable 

case where d =1.40 and d is held at 1.00. The 
xy yx 

additional column on this table represents which force term 

or is increasing without bound while the other term 

is decreasing. It should be noted that some variables were 

given initial values of 1.001 to get comparative plots over 

time, i.e. using 1.25 as an initial parameter in these cases 

caused either M or M to increase or decrease severely 

X y 

immediately . 

Graphically three general forms of instability were 
demonstrated where only one dependent variable was perturbed. 
These are presented in Figures 14, 15, and 16 (trials 1 
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TABLE V 



INITIAL CONDITIONS WITH 

d AT l.Q AND d AT 0 . 1 (STABLE) 
xy yx 



Trial 

# 


Time 


Increment 




1 


0 


to 


50 


0.05 


1.25 


2 


Q 


to 


50 


0.05 


1.00 


3 


0 


to 


50 


0.05 


1.00 


4 


0 


to 


50 


0.05 


1.00 


5 


0 


to 


50 


0.05 


1.25 


6 


0 


to 


50 


0.05 


1.25 


7 


0 


to 


50 


0.05 


1.25 


8 


0 


to 


50 


0.05 


1.00 


9 


0 


to 


50 


0.05 


1.00 


10 


0 


to 


50 


0.05 


1.00 


11 


0 


to 


50 


0.05 


1.25 


12 


0 


to 


50 


0.05 


1.25 


13 


0 


to 


50 


0.05 


1.25 


14 


0 


to 


50 


0.05 


1.00 





I 


_Y, 


Time for 
90% Damp 
Out (Sec) 


1.00 


1.00 


1.00 


2.0 


1.25 


1.00 


1.00 


21.0 


1.00 


1.25 


l.QQ 


21.0 


1.00 


1.00 


1.25 


19.0 


1.25 


1.00 


1.00 


21.0 


1.00 


1.25 


l.QO 


20.0 


1.00 


1.00 


1.25 


18.0 


1.25 


1.00 


1.25 


4.0 


1.00 


1.25 


1.25 


34.0 


1.25 


1.25 


1.00 


6.0 


1.25 


1.25 


1.00 


8.0 


1.25 


1.00 


1.25 


6.5 


1.00 


1.25 


1.25 


33.0 


1.25 


1.25 


1.25 


12.5 
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Figure 13. Single Perturbation, stable phase plot (I at 1.25). 
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TABLE VI 



Time 
0 to 50 

0 to 50 

0 to 50 

0 to 50 

0 to 50 

0 to 50 

0 to 50 

0 to 50 

0 to 50 

0 to 50 

0 to 50 

0 to 50 

0 to 50 

0 to 50 



INITIAL CONDITIONS WITH 

AT 1.4 AND d AT 1.0 (UNSTABLE) 
xy yx 



Increment 


I 

X 




I 


M 


Increasing 
Force Term 


.05 


1.25 


1.00 


1.00 


l.OQ 




.05 


1.00 


1.25 


1.00 


1.00 


X 


.05 


l.QQ 


l.OQ 


1.001 


1.00 




.05 


1.00 


l.QQ 


1.00 


1.001 




.05 


1.25 


1.25 


1.00 


1.00 


^x 


.05 


1.001 


1.00 


1.001 


1.00 




.05 


1.001 


l.QQ 


1.00 


1.001 


M 

y 


.05 


1.00 


1.00 


1.001 


1.001 


^y 


.05 


1.00 


1.001 


1.00 


1.001 


xM 

X 


.05 


1.00 


1.01 


1.01 


1.00 


M 

y 


.05 


1.25 


1.25 


1.25 


1.00 


M 

X 


.05 


1.25 


1.25 


1.00 


1.25 


M 

X 


.05 


1.001 


1.000 


1.001 


1.001 


M 

y 


.05 


1.00 


1.001 


1.001 


1.001 


M 

y 
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through 4)_. The corresponding phase plots are given in 
Figures 17, 18, 19, and 20. Specific values for double 
and triple perturbation did not demonstrate any unusual 
graphic forms of instability and generally followed those 
graphic representations displayed for single perturbations. 

E. INSTABILITY WITH d^^ AT 1.4Q 

Table VII gives the initial values for trial runs holding 

d at 1.00 and d at 1.40. Similar graphic forms were 
xy yx 3 

produced in this investigation as those produced in the 
preceding investigation of instability. Two exceptions are 
given in Figures 21 and 22 (trial 7) and Figures 23 and 24 
(trial 10) . It also should be noted that in this case, 
perturbing M and I together causes M to increase. 

In the previous case perturbing these variables together 
drove upward, an interesting result. 

F . GENERAL OBSERVATIONS 

The cases investigated in sections B and C of this 

chapter clearly indicate that the Dynamic Model of Modern 

Military Conflict does exhibit stability as predicted in 

Chapter 3. Furthermore, it also may behave in an unstable 

manner as indicated in sections D and E above. A closer 

look at the unstable cases also reveals bifurcation relative 

to the values for d and d , initial values for the 

xy yx' 

dependent variables and the resulting increasing term. 
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Figure 17. Single Perturbation, unstable phase plot (I at 1.25). 
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Figure 18. Single Perturbation, unstable phase plot (M at 1.25). 
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Figure 19. Single Perturbation, unstable phase plot (1^ at 1.001). 
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Figure 20. Single Perturbation, unstable phase plot (M at l.OQl). 



TABLE VII 



INITIAL CONDITIONS WITH 









d 

xy 


AT 1.0 


AND d 

yx 


AT 1.4 


(UNSTABLE) 




Trial 

# 


Time 


Increment 


I 

X 


M 

X 


I 

_y 


M 

_z 


Increasing 
Force Term 


1 


0 


to 


50 


.05 


1.25 


1. 00 


1.00 


1.00 


M 

X 


2 


0 


to 


50 


.05 


1.00 


1.25 


1.00 


1.00 


M 

X 


3 


0 


to 


50 


.05 


l.QQ 


1.00 


1.25 


1.00 


M 

y 


4 


0 


to 


50 


.05 


l.OQ 


1.00 


1.00 


1.25 


M 

y 


5 


0 


to 


50 


.05 


1.25 


1.25 


1.00 


1.00 


M 

X 


6 


0 


to 


50 


.05 


1.01 


1.00 


1.01 


1.00 


M 

y 


7 


0 


to 


50 


.05 


1.01 


1.00 


1.00 


1.01 


M 

y 


8 


0 


to 


50 


.05 


1-00 


1.00 


1.01 


1.01 


M 

y 


9 


0 


to 


50 


.05 


1.00 


1.01 


1.00 


1.01 


M 


10 


0 


to 


50 


.05 


1.00 


1.01 


1.01 


1.00 


M 

X 


11 


0 


to 


50 


.05 


1.25 


1.25 


1.25 


1.00 


M 

X 


12 


0 


to 


50 


.05 


1.25 


1.25 


1.00 


1.25 


M 

X 


13 


0 


to 


50 


.05 


1.001 


1.00 


1.001 


1.001 


M 

y 


14 


0 


to 


50 


.05 


l.OQ 


1.001 


1.001 


1.001 


M 
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Figure 21. Double Perturbation, unstable (I and M at 1.01) . 
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Figure 22. Double Perturbation, unstable phase plot (I and 
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Figure 24. Double Perturbation, unstable phase plot (M and I at 1.01) . 



A simple display of this bifurcation is seen when only 

one element (dependent variable), is perturbed as in trials 

1 through 4 in Table VI and VII . Relative to Table VI , d^^ 

is set at 1.40 and d at l.OQ. This gives Y an 

advantage over X in c3 effectiveness, however, it can be 

seen that in trials 1 and 2 of this table, X can overwhelm 

Y (M^ increasing, decreasing) if given an advantage in 

either information (I ) or forces {.M ) . Similarly, relative 

to Table VII with d set at 1.00 and d at 1.40, Y 

xy yx ' 

can be seen to overwhelm X in trials 3 and 4. Evidences of 
bifurcation are a bit more complicated in the remaining trials 
of these tables. 

In Table VI, trials 5, 9, 11, and 12 exhibit this 
bifurcation where y has an advantage in but is over- 
whelmed by X. Similarly in Table VII, trials 6, 7, 8, 13 
and 14 show X being overwhelmed by Y, even though it has 
a C advantage. Trail 10 is not a bifurcation example, but 
as mentioned in section E, does produce an interesting result. 
Specifically, from Table VI, a slight advantage in M 

3 

-and I^ results in Y overwhelming X when Y has a C 

advantage. Whereas from Table VII, the slight advantage in 

M and I results in X overwhelming Y when X has the 

X y 

advantage. This seems consistent with common sense. 
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V. EQUILIBRIUM ANALYSIS 



All the preceding analysis of stability was predicated 
upon a constant equilibrium vector of 1 for ease of calcula- 
tion. Another question that should be answered at this point 
is "what happens to equilibrium values if, for instance, when 
replenish rates (Q^,R^,Q^,Ry) are kept constant and d^^ and/or 
d^^ is allowed to vary?" To answer this question the equi- 
librium equations and initial analysis conditions are brought 
back into the analysis. 



A. METHOD OF ANALYSIS 

First, the original equilibrium equations are set to 0 and 
solved for the replenishment rates. This result is then used 
to solve the same set of equations for equilibrium as a func- 
tion of d and d 

xy yx 



B. EQUILIBRIUM RELATION TO d AND d 

xy yx 

Each equilibrium equation is set to 0 and solved for 

replenishment rates. 

CM -a I -yM -A+Q=0 
XX xe xy ye ' xy ye x x 



0(1) - (1) (1) - 0(1) - 0.5 + Q^ =0 

Q =1.5 

X 



(Eqn 5.1) 



-dl -(6+d)I -(S +b)M -B+R=0 

XX xe xy xy' ye ^ ^xy xy' ye x x 



0(l)-(0 + d )l -(1 + 0) 1-0. 5 + R =0 

X 



^x = ^xy + 1-5 



(Eqn 5.2) 
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0 



-a I -yM +C M -A+Q = 
yx xe 'yx xe yy ye y y 

Q (1) - 1 (1). + Q (1). - 0.5 +Q^ = 0 

Q =1.5 

y 



(Eqn 5.3) 



-(6 +d)I -(3 +b)M -dl -B+R=0 

yx yx xe yx yx xe YY ye y y 

- (0 + d )1 -{1 + 0)1 - 0 (1) - 0.5 + R = 0 

yx y 



R^ = d^^^ + 1.5 (Eqn 5.4) 

y 

utilizing these results the same equations are solved 

again for the equilibrium points starting with some arbitrary 

initial values for d and d / d' and d' respectively 

xy yx ' xy yx ^ ^ 

to fix the replenishment rates. 

CM -a I -yM -A+Q=0 

XX xe xy ye xy ye x x 

0(M )-(l)(I )-0(M )-0.5+Q = 0 
xe ye ye x 



I = Q - 0.5 
ye X 

Using the results of equation 5.1 



XX xe 



xe 



0) 

>1 

H 


= 1 




(6 


+ d ) I 


- (S 


xy 


xy ye 


(0 + 


d ) I 


(6 


xy' ye 


xy 


M 


= -d I 


- 0. 


ye 


xy ye 



(Eqn 5.5) 



xy xy ye 

+ 0) M - ( 
ye 



X 



X 



Using the result of equations 5.2 and 5.5 and substituting 



+ 1.5 for 



M = _d - 0.5 + (d ' + 1.5) 
ye xy xy 



M ^ = 1.0 - (d - d ') 
ye xy xy 



(Ean 5.6) 
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0 



-a I -Y M +CM -A+Q = 
yx xe 'yx xe yy ye y ^y 



- (0>lxe - + <“>V ■ *°y 

' °y - 

Using the results of equation 5.3 



M =1 
xe 



= 0 



(Eqn 5.7) 



-(5 


+ d )I 


(8 


+ b 


)M 


- d I 


- B + R = 0 


yx 


yx xe 


yx 


yx 


xe 


0) 

>1 

>1 

>1 


y y 


- (0 


+ d ) I 
yx' xe 


( 1 


^ “ >-“xe 


- (0) I ^ 

ye 


- 0 . 5 + R^ = 0 



Using the results of equation 5.7 

1.5 - R 



xe 



-d 



yx 



Using the results of equation 5.4 and substituting 

d ' + 1.5 for R 
yx y 

d' + 1.5 - 1.5 

I = ^ 



xe 



yx 



Adding d-d to the numerator 
^ yx yx 

(d - d' ) 

I = 1 - yx yx 



xe 



yx 



C. ASSIMILATION 

It suffices to say then that I and M are not 

ye xe 

effected by changes in d^^ or ^yx* other hand 

varies in the opposite direction to deviations from 
its nominal value. I varies in the opposite direction to 
deviations from its nominal value normalized by ^xy* 
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More specifically, equation 5.6 indicates that force Y 
would require fewer general purpose resources to hold force 
X to 1 provided, force Y increased its enhancement 
beyond its initial value. In a similar fashion equation 5.8 
indicates force X would not require as much information to 

3 

hold force Y to 1 provided force X increased its C 
enhancement beyond its intial value. It should be noted 
that these equations may not be true in general , but rather 
may only be true for the specific initial values used in 
this analysis. 
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VI. SUMMARY AND CONCLUSIONS 



This thesis has described the initial investigation of 
the Dynamic Model of Modern Military Conflict. Following a 
brief description of the model, investigation of stability 
was undertaken using a method of analysis previously prescribed 
for the study of deterministic systems such as ecosystems. 
Through this investigation it has been demonstrated that the 
model does in fact behave as one might expect, given a set of 
somewhat reduced initial parameters. In this respect, this 
research endeavor did not attempt to fully prove or disprove 
the model but rather its purpose was to permit observation of 
the model behavior as an initial step toward further study. 

A. GENERAL OBSERVATIONS 

It can be seen that the Dynamic Model of Modern Military 

Conflict, given a set of asymmetrical initial parameters 

exhibits behavior that corresponds to a status quo (stability) 

or a dominance of one force, combat force if you will, over 

another (instability). . The significance of this observation 

lies in the plausible use of this model in predicting outcomes 

of military conflict. Particularly, the incorporation of 

3 3 

terms relating ro counter C and C enhancements along with 
more traditional terms such as loss rates, replenishment 
rates and Lanchester coefficients, adds a contemporary flavor 
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to this model. The actual employment of this model in 
playing "what if?" type wargames however cannot be undertaken 
as yet without additional research including the remaining 
terms of this model. 

B. IMPLICATIONS FOR FUTURE RESEARCH 

Specifically, the following areas for additional research 
are suggested based on the initial investigation described in 
this thesis; 

1. Location of other equilibrium points besides the 
vector 1 utilized herein. 

2. Determination of an appropriate Lyapunov function 
suitable for predicting global stability. 

3. Prediction of bifurcation through an analytical 
approach. 

4. Further investigation of model behavior relative to 
variations of other initial parameters (i.e. counter C3 
coefficients ) . 



65 



LIST OF REFERENCES 



1. Cushman, John H., "Exercise, Field Test and 
Experimentation and Battle Simulation Approaches," 
Proceedings for Quantitative Assessments of the 
Utility of Command and Control Systems , MITI^ 
Corporation, pp. 37-48, January 1980. 

2. Moose, Paul H. , "Discussion and Notes Regarding 
the Dynamic Model of Modern Military Conflict," 
unpublished notes, p. 1-14V, November-December 1981. 

3. Ibid, p. l.-V. 

4. May, Robert M. , Stability and Complexity in Model 
Ecosystems , Princeton University Press, 1974 . 

5 . Eigenvalues and (optionally) Eigenvectors of a Real 
General Matrix in Full Storage Mode, Copyright by 
IMSL, Inc., 1978. 

6. Hilleary, R. R. , Interactive Ordinary Differential 
Equations Package, W.R. Church Computer Center, User ' s 
Guide to VM/CMS at NPS , April 1981. 



66 



INITIAL DISTRIBUTION LIST 



No. Copies 

1. Defense Technical Information Center 2 

Cameron Station 

Alexandria, Virginia 22314 

2. Library, Code Q142 2 

Naval Postgraduate School 

Monterey, California 9394Q 

3. Department Chairman, Code 74 1 

C3 Academic Group 

Naval Postgraduate School 
Monterey, California 9394Q 

4. Superintendent, Code 39 3 

C3 Curricular Officer 

Naval Postgraduate School 
Monterey, California 9394Q 

5. Captain Brian Van Ormann 1 

AFIT/CIRF 

Wright-Patterson AFB , Ohio 45433 

6. Air University Library 1 

AUL/LD 

Maxwell AFB, Alabama 36112 

7. Professor Paul H. Moose, Code 62Me 5 

Department of Electrical Engineering 

Naval Postgraduate School 
Monterey, California 9394Q 

8. Professor Glenn F. Lindsay, Code 55Ls 1 

Department of Operations Research 

Naval Postgraduate School 
Monterey, California 9394Q 

9. Major Timothy W. Fox, USAF 2 

P.O. Box 35Q0 

MacDill AFB, Florida 33608 



67 



'I 07c Q 9 

Thesis 1 C u 

F6647 Fox 

c,l An investigation of 

the dynamic model of 
modern military conflict 

$6AUC-8» 2 9Q59 



Thesis 
F6647 
c. 1 



1S7G82 

Fox 

An investigation of 
the dynamic model of 
modern military conflict. 



thesF6647 

An investigation of the dynamic model of 




3 2768 000 99861 1 

DUDLEY KNOX LIBRARY 




